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The present application claims priority to, and incorporates by reference, U.S. 
Provisional Patent Application Serial No. 60/456,766, which was filed on March 21, 
2003. 

Background of the Invention 
L Field of the Invention 

This disclosure relates generally to gravity and its uses in applications such as 
drilling and logging. 

2. Background 

A. Subsurface Gravity Background 

Subsurface gravity may have began with pendulum measurements made during 
1826, 1828, and 1854 by Airy (1856) who sought to determine the mean density of the 
earth by measuring the interval vertical gradients between the top and bottom of various 
mine shafts. Similar measurements were made by von Sterneck in 1883 and 1885 and by 
Rische from 1871 to 1902. It was not until the development of the portable gravity meter 
in the 1930s that extensive studies of subsurface gravity became possible. Since then, 
gravity measurements have been made in mine shafts to determine the densities of 
adjacent rocks, to determine the mean density of the earth, and to study anomalous 
vertical gradients caused by the large positive density contrasts associated with ore 
bodies. The concept of a borehole gravity meter for gravity logging of wells was first 
proposed in the 1950s, and development of the high-precision borehole gravity meter 
(BHGM) began in the 1960s. 

Subsurface gravity measurements (e.g. in boreholes or in mines) are used to detect 
anomalous density structure or lateral density variations. These density variations are 
typically due primarily to stratigraphic, structural, or diagenetic effects that cause 
subsurface iso-gravity contours of equal gravity to depart from the horizontal. However, 
interpreting density variations is difficult, in part, because density models are not 
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unique — in other words, different models can be used to describe the same density 
variations. In theory, an infinite number of density-volume model configurations can be 
devised to generate the same gravity anomaly. But, in many cases, a uniform, horizontal 
layered earth can be assumed as a useful model because formations surrounding gravity 
5 measurements are typically lateral or nearly-so and possess relatively-uniform densities 
in lateral directions. In such areas, subsurface gravity data are easily converted to 
accurate and unique interval density profiles. 

Lateral density variations may be significant where folded strata, faults, 
unconformities, intrusions, or lateral variations in lithology, porosity or pore fluids (due 
to selective depositional or post-depositional processes) intersect or occur within 
detectable distances of a borehole. Under these circumstances, analysis of borehole 
gravity data is more difficult because equal density surfaces generally are poorly known 
and may be complex in shape. 

While most other geophysical logging tools only sample in the local area of the 
borehole, the classical BHGM log samples a large volume of rock. Consequently, it is 
not significantly influenced by drilling mud, fluid invasion, hole rugosity, or formation 
damage that surrounds all wells to some degree. It is the ability of the BHGM log to 
yield a direct measure of in situ density and porosity and to characterize pore content that 
forms the basis of its application to oil and gas exploration and production. 

Borehole gravity surveys have proven to be worthwhile at least because of their: 
(1) High relative or absolute accuracy; 
25 (2) Direct density response; and 

(3) Ability to investigate a large volume of formation. 

Borehole gravimetry is also useful due to its ability to measure apparent density. 
Traditional logging tools only measure into a formation a few inches, whereas borehole 
30 gravity samples tens to hundreds of meters into the formation. One aspect of the BHGM 
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that makes it an attractive logging tool in the petroleum industry is its ability to detect the 
contacts between gas, oil, and water at large distances from the borehole. It also can do 
this through multiple casing strings and formation damage, which decreases the 
effectiveness of the competing pulsed neutron tools for fluid saturation monitoring, as 
does low salinity water. 



B. Physics of Subsurface Gravity 
Gravity exploration utilizes Newton's Law of Universal Gravitation, 



Force = G 1 . (1) 
r 



This law states that, between any two massive objects, there is a force, which is directly 
proportional to the product of their masses and inversely proportional to the square of the 
distance between them. Thus, the force of attraction is larger for larger masses but 
decreases rapidly for increasing distances between the masses. The constant of 
proportionality, G, is the Universal Gravitational Constant. A gravity sensor (or gravity 
meter) uses a very sensitive mass balance mechanism to measure the force of gravity 
acting on a test mass (often called a proof mass). 



Using Newton's Second Law of Motion, the force on a test mass m 2 can be 
expressed as 

Force = m 2 g , (2) 

where g is the gravity or gravitational acceleration experienced by the test mass due to 
the force of gravity. Equations (1) and (2) can be combined to give: 

m 2 r 

Thus gravity is expressed in units of force per unit area or acceleration. In 
geophysics, the units of centimeters/second is called "gal ", after Galileo. Land, airborne 
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and marine gravity surveys are typically recorded in units of milli-gal or mgal and 
borehole gravity surveys are typically given in units of micro-gal or /ugal. 

Although gravity is a tensor, most exploration gravity meters are only sensitive to 
the vertical component of gravity, g z , which is in the direction of the local plumbline. 
Therefore, the acceleration experienced by a gravity meter test mass is 



r 2 



g 2 =G^cos<p, (4) 

where q> is shown in FIG. 1 . 



Usually gravity exploration involves the spatial distribution of rocks whose 
densities and volumes are partly known or can be inferred. Because mass = density x 
volume, we can rewrite equation (4) as 

g z =G^cos<p , (5) 
r 

where p is the density and Vis the volume. For practical purposes, we integrate Equation 
(5) over the entire volume, 



g 2 =G\£;Cos<pdV. (6) 

Equation (6) calculates the vertical component of gravity at any point due to any mass, 
where mass is defined in terms of density and volume. Therefore density models can be 
constructed using Equation (6) to simulate proposed geologic models and their gravity 
fields so that the models can be compared with actual measured gravity values. There are 
many commercially-available software packages available for doing this. Many of them 
are based on an algorithm developed by Talwani (1965). 
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If the earth is assumed to be non-rotating and perfectly spherical with a radially 
symmetric density distribution, with a radius of R and a mass of M, then, from Equation 
(3), the gravity on the earth's surface will be 
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and the free-air gradient of gravity at the surface of the earth will be 

* -|«<% (8) 
or r = R 3 

where p R is the bulk mean density of the earth. At any point inside the earth at a radius 
of r from the earth's center, the mass of a sphere of radius r is given by 

m r =47rj r Q p(r)r 2 dr (9) 

where p r is the internal density as a function of r. Substituting m r in Equation (9) for M 
in Equation (7) gives the gravity at this point (Benfield, 1937) 

g=— j 0 P(r)r 2 dr (10) 

because the net attraction of a spherical shell between r and is zero (Ramsey, 1940). 
The vertical gradient of gravity at this point is 

which reduces to 



%L = 4nGp-\nGp ri (12) 

30 where p is the density of an infinitesimally thin spherical shell of radius r and p r is the 
mean density of the interior sphere of radius r. See Airy (1856), Miller and Innes (1953), 

25397625.1 

-6- 



Gutenberg (1959), Hammer (1963), and Beyer (1971) for further discussions. From 
Equation (8) we see that the second term in Equation (12) is the free-air vertical gradient 
of gravity for a non-rotating spherical earth. 



In order to predict the gravitational field of the earth precisely at any point on the 
earth, we must know and correct for the shape and density distribution with the greatest 
possible accuracy. The earth is not actually a sphere because it rotates and thus bulges at 
the equator and flattens at the poles. Its shape can be closely approximated by an oblate 
spheroid with an eccentricity of 1/297. The rotation and general ellipsoidal shape of the 
earth can be taken into account by replacing the second term in Equation (12) with the 
normal free-air vertical gradient, which can be determined by using the spheroid model 
based on a best-fit reference model of gravity at mean sea level as a function of latitude 
This reference standard model is established by the International Union of Geodesy 
and Geophysics. The accepted model was last updated in 1967 and is 



This reference standard model says that gravity varies at mean sea level from 
978,030.90 mgal at the equator to 983,884.10 gal at the poles, which is a total range of 
5,853.20 mgal. Note that the gravity is larger at the poles than at the equator because the 
equator is farther from the earth's center of mass than are the poles. The normal free-air 
vertical gradient is 



£i967 =[9.7803090 + 0.058552sin 2 ^-5.70xl0" 5 sin 2 2^] ms~ 2 



(13) 



or in ft s 



g l967 =[32.0875312 + 0.192099sin 2 ^-1.87xl0- 4 sin 2 2^] fts~\ 



(14) 




0.0941 12 - 0.0001 34sin 2 0-O.134xlO~ 7 h 



(15) 
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where h is the elevation in feet. The normal free-air gradient of gravity varies from the 
equator to the poles by less than 0.2% and with elevation by about 0.01% per 1,000 feet 
or 0.05% per kilometer (see Hammer, 1970). These variations are very small for gravity 
surveys made in boreholes and can be approximated by 

5 

F = |£ = 0.09406 mGal/ft . (16) 

When this is substituted for the second term in Equation (12), we obtain the 
vertical density p of the laterally adjacent rocks for the case in which the earth possesses 
10 a radially symmetrical distribution of density 

— = AnGp -F = AnGp - 0.09406 ( 1 7) 

dr 



15 



or, changing from elevation, h, to depth, z, we have 

Ag = FAz-AxGpAz (18) 

or 

Ag = 0.09406 Az-AnGphz (19) 

20 Equation (18) is a fundamental equation of borehole gravity. The4;rGp term 

applies to an infinitely extended horizontal layer of thickness Az . According to Equation 
(18), gravity increases downward at a rate determined by the difference between the free- 
air vertical gradient F, which is essentially a constant, and a gradient of opposite 
sign, AnGp, that varies as the density of the adjacent rocks change. The positive term is 

25 due to an increase in gravity downward caused by closer approach to the center of the 
earth and the negative term is twice the attraction of an infinitely extended horizontal 
layer of thickness Az . Therefore, increases in layer density correspond to decreases in 
the interval vertical gradient Ag/Az, and vise versa. Thus it is often possible to 
accurately determine very small variations in the density of rocks bracketed by different 
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Az intervals with little or no analysis beyond the simple reduction of the basic gravity 
and depth measurements, which is one of the great strengths of borehole gravity. 

C. Present Borehole Gravity Logging Technology 
5 FIG. 2 illustrates a conventional borehole gravity survey. At present, borehole 

gravity logging is conducted using a BHGM, which consists of a single gravity sensor 
located in a logging sonde (10). The BHGM logging sonde is lowered down the well (12) 
on the end of a wireline (14). All the measurements are taken at stationary locations 
because any gravity sensor acceleration will be measured in addition to the gravity (due 
10 to the Principle of Equivalence). Measurements are taken at (16) and (18) in order to 
determine the apparent density of the intervening layer. The first measurement is taken at 
(16). The gravity measured at (16) is due to the mass of the intervening slab of density 
p and of thickness Az , as given by Equation (18), plus the gravity due to the mass of the 
rest of the universe, 

15 

g m = InGptsz + g{universe) . (20) 

The second measurement is taken at (18). Similarly, the gravity measured at (18) is due 
to the mass of the intervening slab of density p and of thickness Az plus the gravity due 
20 to the mass of the rest of the universe, 

£(i2) = -ZnGpAz + g(universe) (21) 
The difference between these two measurements is 

25 

A S = S(i2)-£(u» (22) 
Ag = -4xGpAz, (23) 

30 which gives 
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= -1 Ag 
AkG Az 



Therefore, the density of a layer of rock intersected by a well can be directly determined 
by measuring the gravity on either side of that layer and the distance between those 
5 gravity measurement points. 

Equations (18) and (24) are not strictly correct when applied to the real earth. 
Departures of the earth's surface from an ellipsoid and lateral density variations in the 
subsurface contribute to the variation of gravity in the subsurface. Unwanted or 
10 extraneous accelerations caused by topography or mass disturbances connected with the 
well are usually negligibly small and can be ignored or corrected-for with sufficient 
accuracy. The corrections for terrain can be found in most general geophysical textbooks 
{e.g. Dobrin, 1976; Telford et al., 1976) and will not be discussed here. 

15 Present commercially-available borehole gravity meter instruments are only 

capable of making measurements to a deviation of less than 14 degrees from vertical. 
Some prototype borehole gravity meters have been developed that can make 
measurements up to a deviations of as much as 115 degrees from vertical. FIG. 3 is an 
illustration of such a borehole gravity meter (20). In principle, the gravity sensor (22) is 

20 of any design that can make gravity measurements sufficient to meet its applications and 
that can be placed in a small enough housing to accommodate its use in a borehole 
logging environment - see for example Ander et al. 9 1999b. In practice, such sensors 
primarily consist of either a metal or a fused quartz relative spring sensor element. The 
gravity sensor (22) housing is gimbaled about its horizontal axis (24) and its vertical axis 

25 (26). The gravity sensor (22) is oriented to vertical by activating motors (30A) and (30B) 
that are connected to the sensor axes (24) and (26) through a series of gears (28A), (28B), 
(28C) and (28D). The motors are controlled by commercially-available high precision 
levels usually located in the gravity sensor housing. 
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Once borehole gravity instruments are capable of routinely making measurements 
at significant deviations, then corrections for well deviation can become important. 
Corrections for borehole deviation is generally straight forward. The correction is 

Az^=Az„,cos^ (25) 

It converts the measured depth to true vertical depth, where 6 is the deviation from 
vertical. Rivero (1971) and Fitchard (1981) have developed more sophisticated 
corrections for the case where the well deviation involves doglegs. If corrections are not 
made for deviations, then the interval gravity gradient can be erroneously low, calculated 
interval densities can be erroneously high, and interval porosities can be erroneously low. 

If a well is significantly deviated, then the theoretical latitude-dependent 
horizontal gradient of total gravity and any anomalous horizontal gradients of total 
gravity will contribute to the gravity variations measured down the well. The theoretical 
latitude-dependent gradient does not exceed 0.2 figal/ft and applies only to the north- 
south component of deviation. The actual value is given by 0.8 122 sin 2^ mgal/km, 
where (j> '\s latitude (Nettleton, 1976, p. 80-81). Anomalous horizontal gradients in total 
gravity occasionally are greater than 1.9 to 3.8 jugal/ft. Values of the theoretical latitude- 
dependent horizontal gravity gradient together with estimates of anomalous horizontal 
gradients taken from surface gravity maps can be analyzed with hole azimuth and hole 
angle data from the well directional survey to determine if corrections for horizontal 
gravity gradients are necessary. Using present borehole gravity meter technology, in the 
vast majority of cases, corrections for these effects are unnecessary or are very small. 

Once borehole gravity measurements can be made in horizontal or near horizontal 
wells, then the use of Equation (24) breaks down because Az -» 0 and there is no longer a 
significant vertical gradient between successive gravity measurements. In the horizontal 
well logging case, a gravity logging survey can be treated like a horizontal land gravity 
profile. 
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D. Shortcomings 
Despite the technology outlined above, significant shortcomings remain relating 
to gravity measurement and the use of gravity measurements in applications such as 
logging and drilling. Namely, conventional technology does not provide for the ability to 
5 (a) perform gravity well logging using arrays of gravity sensors; (2) create density 
pseudosections; (3) make gravity measurements while drilling; and (4) steer a drill bit or 
other apparatus using gravity measurements. Techniques of this disclosure, however, 
address these shortcomings, as discussed in detail below. 

10 Summary of the Invention 

The technology described in this disclosure can be applied to a vast array of 
applications, including drilling and logging applications. 

More particularly, techniques are discussed that allow one to make gravity 
15 measurements in wells, both during the well-drilling phase and afterward. Gravity 
measurements can be made during well logging. Gravity measurements can be made 
using one or more gravity sensors attached to a drill string or a logging string. The 
general process of making sensor measurements during drilling is sometimes referred to 
as Logging-While-Drilling (LDW) or Measurement- While-Drilling (MWD), and 
20 embodiments of the present disclosure are referred to as Gravity-Measurement- While- 
Drilling (GMWD). 

Embodiments of this disclosure also involve the use of GMWD to steer a drill bit 
and the use of subsurface iso-gravity contour maps. Still other embodiments involve 
25 making multiple gravity measurements using multiple gravity sensors attached along a 
logging string or otherwise forming an array. In different embodiments, such gravity 
array measurements can be made simultaneously to provide a significant amount of 
useful data in a short period of time. 
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In one respect, the invention involves an apparatus including a plurality of gravity 
sensors spaced a known distance apart to form a gravity sensor array, the array being 
coupled to a subsurface logging apparatus. The logging apparatus can include a logging 
sonde, coil tubing, or wireline. Sensor drift associated with one gravity sensor can be 
5 time-correlated with a sensor drift of another gravity sensor. The sensor drifts of all the 
gravity sensors can be time-correlated with one another. The sensor drift of one or more 
of the gravity sensors can be corrected using sensor hardware or software. Different 
gravity sensors can be spaced at different distances from one another. The array can be a 
linear, vertical array. 

10 

In another respect, the invention involves a method of logging using multiple 
gravity sensors. Gravity sensors are spaced a known distance apart to form a gravity 
sensor array. The gravity sensor array is placed into a subsurface, and the gravity sensor 
array is used to make gravity measurements of the subsurface. The gravity sensors can 
1 5 make gravity measurements simultaneously. 

In another respect, the invention involves a method in which an apparent density 
is calculated. The apparent density is plotted as a function of position and gravity point 
separation to form a density pseudosection. The density pseudosection indicates how the 
20 apparent density varies with location and depth. The method can also include forming 
contour lines associated with the apparent density. 

In another respect, the invention involves an apparatus including a gravity sensor 
coupled to a drilling device, the sensor being configured to make a subsurface gravity 

25 measurement after the drilling device begins drilling a hole but before the hole is 
completed at a time during which the drilling device is temporarily halted. The gravity 
sensor can be coupled to an instrument package. The instrument package can be coupled 
to a drilling pipe, which is coupled to a drill bit. The sensor can be configured to make a 
subsurface gravity measurement when a drilling pipe is added to a drill string. The 

30 gravity sensor can be spaced a known distance from one or more other gravity sensors to 
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form a gravity sensor array. The gravity sensor can be a metal zero-length spring gravity 
sensor. The gravity sensor can include a clamping mechanism configured to withstand 
vibrations due to drilling. 

5 In another respect, the invention involves a method involving coupling a gravity 

sensor to a drilling device. A hole is drilled using the drilling device. Drilling is 
temporarily halted. A gravity measurement is made while drilling is temporarily halted, 
and drilling of the hole is then resumed. The method can also include using the gravity 
measurement to steer the drilling device. The method can also include using the gravity 

10 measurement to steer the drilling device into and out of salt structures. The method can 
also include using the gravity measurement to steer the drilling device into porous 
structures. The method can also include using the gravity measurement to steer the 
drilling device into compartments formed by up-thrown or down-thrown faulted blocks. 
The method can also include using the gravity measurement to steer the drilling device 

15 into or around man-made structures. The method can also include using the gravity 
measurement to verify a course of the drilling device. The hole can be a well. The 
method can also include using the gravity measurement to locate porosity. The method 
can also include using the gravity measurement to profile a salt interface. The method can 
also include using the gravity measurement to log carbonate units. The method can also 

20 include using the gravity measurement to detect fractures. The method can also include 
using the gravity measurement to determine reef proximity. The method can also include 
using the gravity measurement to detect fluid contacts. The method can also include 
using the gravity measurement to delineate lateral facies changes. The method can also 
include using the gravity measurement to determine a drill hole position. The drill hole 

25 position can be an optimized drill hole position relative to a contact. 

In another respect, the invention involves a method comprising steering a drill bit 
using gravity measurements made while drilling. The method can also include 
comparing the gravity measurements with modeled gravity. The modeled gravity can 
30 include a subsurface gravity map. 
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In another respect, the invention involves a method for steering a drill bit. A 
geological model is provided of the ground. A gravity calculation is performed based on 
the geological model. A subsurface gravity map of the ground is created based on the 
5 gravity calculation, and the subsurface gravity map is used to steer a drill bit. Using the 
subsurface gravity map to steer a drill bit can include: (a) making a gravity measurement 
during a temporary halt of drilling; (b) using the gravity measurement and the subsurface 
gravity map to define a drilling course; and (c) sending the drilling course to a drilling 
guidance system. The method can also include using a geological model to define the 
10 drilling course. 

Other advantages and features will be apparent to those having ordinary skill in 
the art through reference to the description below. 

15 Brief Description of the Drawings 

The following drawings form part of the present specification and are included to 
further demonstrate certain aspects of this disclosure. The invention may be better 
understood by reference to one or more of these drawings in combination with the 
detailed description of specific embodiments presented here. Use of identical element 

20 numbers among the drawings denotes like or identical structures. The drawings are not 
to scale and are included for illustration only. 

FIG. 1 shows a gravitational attraction of a point mass. 

25 FIG. 2 shows conventional borehole gravity logging. 

FIG. 3 shows a gravity sensor for making gravity measurements according to 
embodiments of this disclosure. 
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FIG. 4 shows gravity well logging using an array of multiple gravity sensors in a 
logging sonde, in accordance with embodiments of this disclosure. 

FIG, 5 shows an example density pseudosection method according to 
5 embodiments of this disclosure. 

FIG. 6 shows an example density pseudosection in accordance with embodiments 
of this disclosure. 

10 FIG. 7 shows gravity measurement while drilling in accordance with 

embodiments of this disclosure. 

Description of Illustrative Embodiments 
The reader is notified that the techniques described in this paper are examples 

15 only. The description is drawn to particular embodiments. Those embodiments, 
however, should not be read as limiting the scope of the invention. The invention is 
defined legally by the claims that issue. For example, claims may not include all the 
features described in conjunction with an embodiment; in that case, the claim is broader 
than the embodiment. Likewise, claims may include different combinations from 

20 different embodiments. Those having ordinary skill in the art will recognize that changes 
can be made to the embodiments listed here without departing from the spirit and scope 
of the disclosure and the spirit, scope, and legal coverage of the claims. 

Different embodiments of this disclosure involve the following concepts: 
25 (1) gravity well logging using gravity sensors arrays; (2) density pseudosections; (3) 
Gravity Measurement While Drilling (GMWD) using single or multiple gravity sensors; 
and (4) geosteering using GMWD. 
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1. Gravity Well Logging Using Multiple Gravity Sensors 

The quality of borehole gravity logging surveys depends on reducing the errors in 

Ag and Az measurements small enough to ensure that calculated interval densities have 

the required precision. 

5 

According to Caton (1981), the effect of Ag error on p is given by 



10 



4? 



Az 



(26) 



and the effect of Az error on p is given by 



Perror 



Az A Az 



(27) 



1 5 The probable error in apparent density, P for borehole gravity measurements is 



20 



_ 0.9539^h 2 Sg 2 +Ag 2 Sh 2 
AnGh 2 



(28) 



where Sg is the uncertainty in the gravity measurement and 8h is the uncertainty in the 
height measurement respectively. 



Note that the uncertainty in formation density increases as Az decreases. A 
typical borehole gravity operation acquires data at an accuracy of between 1 to 5 pigal. 
For a measurement uncertainty of 1 jugal, there is an approximate error in formation 

25 density of less than 0.01 g/cm 3 for Az larger than about seven feet. Presently, gravity 
data is seldom collected with Az of less than 10 feet. An important source of density 
uncertainty comes from relative depth inaccuracies. Typical wireline operation acquires 
downhole position at an accuracy of between 0.1 and 0.2 feet. The formation density 
inaccuracy approaches 0.02 g/cm 3 at a Az of about 16 feet. At larger downhole station 

30 spacing, density errors can be under 0.01 g/cm 3 . 

25397625.1 

-17- 



According to embodiments of FIG. 4, an array of multiple gravity sensors (20) are 
placed at intervals in or on a well logging conveyance such as a logging sonde (10), coil 
tubing, wireline, autonomous vehicle, etc. By precisely knowing the intervals between 
5 the multiple gravity sensors (20) along the logging conveyance, then Az error can be made 

negligible compared to Ag error . As a result, p error remains low for much smaller Az . 

There are several components of Ag error . See Ander et al., 1999a for a study of 
the various components of Ag error and their comparative magnitudes. The introduction of 
10 multiple sensors makes it advantageous for the sensor drift for each of the gravity sensors 
to be time-correlated with all the other gravity sensors in order to not increase Ag error . 

Currently, the component of Ag error due to drift is reduced by either correcting for drift or 

by reducing the sensor drift through hardware and/or software improvements. See 
textbooks, for example Dobrin (1976) or Telford et al. (1976), for details of how to 
15 perform a gravity sensor drift corrections. 

Because Az error can be made negligible compared to Ag error , this allows one to 
place gravity sensors close and still maintain acceptable levels of p error . The closest 
spacing between two gravity sensors is limited only by sensor design. 

20 

One example of an array of multiple gravity sensors in a logging device can 
involve 100 gravity sensors, each 30 cm apart, forming a 30-meter-long array of gravity 
sensors down a borehole. With such a system, one can examine apparent densities given 
Az values ranging from 30 cm up to 30 meters in 30 cm increments. The total number of 
25 different apparent densities for such an array is 4,950. 

In other embodiments, a different number of gravity sensors can be used in an 
array. For example, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20, 21, 22, 
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23, 24, 25, 26, 27, 28, 29, 30, 31, 32, 33, 34, 35, 36, 37, 38, 39, 40, 41, 42, 43, 44, 45, 46, 
47, 48, 49, 50, 51, 52, 53, 54, 55, 56, 57, 58, 59, 60, 61, 62, 63, 64, 65, 66, 67, 68, 69, 70, 
71, 72, 73, 74, 75, 76, 77, 78, 79, 80, 81, 82, 83, 84, 85. 86, 87, 88, 89, 90, 91, 92, 93, 94, 
95, 96, 97, 98, 99, 101, 102, 103, 104, 105, etc. can be used. 

5 

Although the spacing between individual sensors can be constant, it need not be. 
Each sensor in the array can make gravity measurements simultaneously, separately, 
according to zones, or in any pattern desired by the user. For example, in one 
embodiment, each sensor in an array may make a measurement at a given time or a given 
10 time range. In a different embodiment, a group of sensors defining a first zone may make 
a measurement followed by a group of sensors in a different zone. With the benefit of 
this disclosure, those having ordinary skill in the art will recognize that the sensor 
measurements can be apportioned according to need and desire. 

15 Although shown as a linear, vertical, array in FIG. 4, gravity sensor arrays can 

take different forms as needed. For example, the array could form curves, or sensors in 
the arrays may be offset from one another. 

With the benefit of this disclosure, those having ordinary skill in the art will 
20 recognize that data from the gravity sensor array can be collected, analyzed, and used in 
various ways. For example, data from the array can be collected and analyzed at a central 
location and conveyed locally or remotely through the use of computer networks such as 
the Internet. 

25 2. Density Pseudosections 

If x is the number total of gravity points in a down-hole set of gravity data point, 
then there is an apparent density p i } value, given by: 

p u = a r T' J U = l,2,-..x (29) 
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for each (i, j) combination of gravity points with separations of Az. . . There are 



(30) 

2 

total number of p (j for x number of gravity points. 

Density Pseudosections, illustrated in FIG. 5 for an array of 7 gravity points, can 
be used to display the apparent density p tJ values associated with each Ag.j/Az. j. 

Density Pseudosections can be compared to electrical pseudosections that are applied to 
electrical geophysics. (Marshall and Madden, 1959). 

A Density Pseudosection is a plot of apparent density as a function of position and 
gravity point separation. A Density Pseudosection indicates how the apparent density 
varies with location down the well and horizontal distance away from the well. The 
apparent density can be converted into models by geophysical inversion (see for example, 
Blakely, 1996). The p. . can be plotted at a point horizontally to the right (or left) of the 

well at the midpoint between the i th and j th gravity points. Values of apparent density, 
p. j , for each (/, j) pair of gravity points can be plotted on a horizontal line at the point of 

intersection of two lines drawn from the base line of the gravity points in the well at an 
angle a , starting at the midpoints between each set of gravity points. 

For convenience, a = 45° in FIG. 5 and FIG. 6. In this way, each apparent 
density value p i} appears at a point midway between the two gravity points, at a 

horizontal distance to the right (or left) which increases with the n value for the spread of 
gravity points. So in FIG. 5, the apparent density, p 15 , which is obtained from gravity 

point #1 and gravity point #5, is plotted 45° to the right of 1 and 5. Apparent densities for 
other pairs of gravity points are plotted in a similar fashion. Each column of apparent 
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density points can be numbered sequentially, e.g. n = 1, 2, 3, 4, 5, 6 The larger the 

value on n the deeper is the apparent density data point into the formation. Note that 

x = n-\. (31) 

5 

Finally, contours can be drawn on the p i } values. FIG. 6 illustrates an example of 

a contoured Density Pseudosection. The result is a form of two-dimensional plot in 
horizontal section. Notice that a Density Pseudosection represents a slice through a 
10 volume of revolution centered on the well. 

3. Gravity Measurement While Drilling 

FIG. 7 illustrates Gravity Measurement While Drilling (GMWD). In general, 
GMWD involves one or more gravity sensors coupled to a device for forming a hole or 
15 opening in a surface. In general, the gravity sensor(s) make subsurface gravity 
measurements after the drilling has begun but before drilling is finished, the 
measurements being made typically at a time during which drilling is temporarily halted. 

In preferred embodiments, and according to FIG. 7, GMWD involves one or more 
20 gravity sensor(s) (20) located in operative relationship with one or more instrument 
packages (32) located in operative relationship with a drilling pipe (34) associated with a 
drill bit (36). In preferred embodiments, the sensors (20) are located in, on, or near the 
instrument packages (32). The instrument packages (32) are preferably in, on, or near the 
drilling pipe (34). The drilling pipe (34) is preferably above the drill bit. 

25 

The instrument package may also contain other devices located above and/or 
below various gravity sensors, such as instruments, other sensors, electronics, thermal 
cooling devices, communications technology, etc., as required. The instrument package 
shown in FIG. 7 is for illustrative purposes only and shows four gravity sensors (20) 
30 separated by compartments containing other devices (38). Multiple sections of drilling 
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pipe containing gravity sensor arrays can be coupled together to produce continuous 
gravity sensor arrays that are much longer than a single length of drilling pipe. 

Conventional technology does not take gravity data while actually drilling a drill 
5 hole (40) through rock (42). One reason is that any sensor accelerations would couple 
with the true gravity measurement and thus could give false readings (the Principle of 
Equivalence of acceleration and gravity). A second reason is that vibrations associated 
with drilling can overpower the gravity sensor and damage or destroy it. 

10 According to embodiments of this disclosure, gravity measurements can be taken 

while a hole is being formed, but at time in which actual drilling is paused. For example, 
drilling can be paused whenever another section of drilling pipe is added to a drill string. 
Drilling pipe typically comes in 30-foot sections. When drill pipe is added, the drilling 
operation is paused for several minutes. During that time, when the drilling is paused, 

15 one or more gravity sensors (e.g., a gravity array if multiple sensors are used) can take 
gravity data. 

The gravity data can be stored electronically and, then transmitted up to the 
surface at any convenient time. The data storage requirements for gravity data from even 
20 a large array of gravity sensors is small and can be transmitted at a low baud rate. Data 
transmission schemes are commercially available and well known in the art, such as mud 
pulsing (acoustic vibrations of the drilling mud). 

If a GMWD device has only one gravity sensor, then the smallest Az would 
25 typically be 30 feet, the typical length of each drilling pipe section. But with a GMWD 
that includes many gravity sensors in each section of drilling pipe, then much shorter Az 
can be obtained using gravity sensor arrays of any desired length. As discussed above, 
by knowing the intervals between multiple gravity sensors (20) along the drilling string, 
then Az error can be made negligible compared to Ag error . As a result, p error remains low 
30 for much smaller Az . Because Az error can be made negligible compared to Ag error , this 
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allows gravity sensors to be placed as close as desired and still maintain acceptable levels 
of p error . The closest spacing between two gravity sensors is limited by sensor design 

only, and arrays can be as long as desired. Density Pseudosections can be applied to the 
apparent density data that is collected with such GMWD gravity sensor arrays. 

5 

As stated earlier, there are several components of Ag error . See Ander et al., 1999a 
for a study of the various components of Ag error and their comparative magnitudes. In 
preferred embodiments, the sensor drift for each of the multiple gravity sensors is time 
correlated with all the other gravity sensors in order to not increase Ag error . Currently, 

10 the component of Ag error due to drift is reduced by either correcting for drift or by 

reducing the sensor drift through hardware and/or software improvements. See 
textbooks, for example Dobrin (1976) or Telford et al. (1976), for details of how to 
perform a gravity sensor drift corrections. 

15 A gravity sensor for GMWD is preferably capable of high precision, robust 

enough to withstand a drilling environment, and possess a small drift rate (although 
software solutions do exist for correcting for larger drift rates). A GMWD gravity sensor 
preferably has a sturdy sensor clamping mechanism in order to withstand vibrations due 
to drilling. Commercially-available gravity sensors have the necessary attributes — for 

20 example, metal zero-length spring gravity sensors are acceptable. Those having ordinary 
skill in the art will recognize that other sensor and associated equipment can be utilized 
for GMWD applications such as those discussed here. 

Applications for GMWD are vast and can include obtaining borehole gravity 
25 measurements in a drilling environment and/or using gravity measurements to aid in geo- 
steering of an apparatus such as a drill bit into a proposed target, either geological, 
structural, or man-made.The following is a non-limiting list of example applications for 
GMWD: 

1 . Locating missed porosity away from the well while drilling. 
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2. Profiling salt interfaces while drilling. 

3. Logging carbonate units while drilling. 

4. Logging sandstone units while drilling. 

5. Detecting fractures while drilling. 

6. Determining reef proximity while drilling on a reef target while drilling. 

7. Looking ahead through salt while drilling. 

8. Detecting the existence and geometry of nearby fluid contacts while 
drilling. 

9. Geosteering: 

a. Into and out of salt structures. 

b. To delineate lateral facies changes. 

c. To optimize drill hole position relative to a fluid, geological, or 
man-made contact (e.g., a caprock, an oil-water contact, a gas- 
water contact, a gas-oil contact). 

d. Into porous structures, e.g. porous carbonate structures, fracture 
zones, sand units, etc. 

e. Into compartments formed by up-thrown or down-thrown faulted 
blocks. 

f. Into or around man-made structures, e.g. tunnels, underground 
facilities, underground structures, etc. 

Other applications will be apparent to those having ordinary skill in the art, with 
the benefit of this disclosure. 

In each of the above applications, the operation of the GMWD apparatus can be 
similar. A drilling apparatus is operated to begin formation of an opening. Gravity 
measurements are taken using one or more gravity sensors coupled to the drilling 
apparatus. In preferred embodiments, the actual measurements are taken at a time in 
which a drill bit has paused, but at a time still during the overall process of forming the 
hole. Gravity measurements are taken using one or more gravity sensor arrays. In 
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preferred embodiments, gravity sensor arrays are utilized, with each gravity sensor being 
one or more known, predetermined distances from one another. Following the gravity 
measurements, drilling continues, and other rounds of gravity measurements can be taken 
as needed or desired. Data resulting from the GMWD can be stored or transmitted as 
required. In preferred embodiments, the data is transferred to the surface for analysis or 
use in, for example, forming Density Pseudosections. In different embodiments, the data 
obtained while drilling can be transmitted to one or more remote users using computer 
networks such as the Internet. Accordingly, data can be analyzed, viewed, or used 
remotely offsite. Likewise, offsite monitoring of the GMWD operations can occur. 

4. Geosteering Using GMWD 

In one embodiment, a drill bit is steered using gravity measurements made while 
drilling. The gravity measurements can be compared with modeled gravity, such as 
gravity modeled in a subsurface gravity map. 

In one embodiment, a GMWD instrument may make measurements of gravity and 
transmit the measurements to the surface via, for example, acoustic vibrations of the 
drilling mud. This data may be combined with drilling location information and then be 
fed into a computer software system designed to analyze such data, called here a 
Comparative Scenario System (CSS). 

In one embodiment, the CSS is simply a computer system configured to analyze 
gravity data using, for example, known algorithms and/or mathematical techniques. In a 
preferred embodiment, however, the CSS is an artificial intelligence computer system 
that may use a combination of gravitational potential inverse theory, particularly gravity 
ideal body analysis, forward modeling, and/or numerical analysis techniques {see, e.g. 
Ander and Huestis, 1987; Huestis and Ander, 1983; Blakely, 1996; Parker, 1975) to 
rapidly evaluate the GMWD gravity data by comparing it to a set of geological scenarios, 
which have been converted into two-dimensional or three-dimensional gravitational field 
subsurface iso-gravity contours. Through this comparison, the CSS may calculate a 
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drilling course correction that it may automatically implement through a drilling down- 
hole guidance system. 

In one embodiment, the geological scenarios and subsurface iso-gravity contours 
5 may be generated by a team of specialists, e.g., geologists, geophysicists, engineers, 
technicians, etc. responsible for the prospect, the reservoir analysis package, or 
intersecting the target. The more known and suspected geological information and/or 
target information the team can place into the models, the greater accuracy the subsurface 
iso-gravity contours will have, and the greater will be the value of the GMWD data to 
1 0 guide the drilling effort. 

Based on the geological/geophysical input, the system may calculate three- 
dimensional gravitational subsurface iso-gravity contours and associated error envelope 
that the drilling unit should stay within for the entire course of the drilling operation. The 
15 system may be designed to give the team real-time feedback by transmitting the data and 
analysis back to the team. 

As long as the GMWD data stays within the error envelope predicted by the 
subsurface iso-gravity contours based on the geological scenarios, the system may 
20 continue automatically. But if the scenarios prove to be substantially incorrect, the 
system may inform the team and ask for a scenario correction. The system permits the 
team to intervene and make remotely-accessed, real-time course corrections to the 
drilling operation based on comparative analysis of the actual data obtained near the drill 
bit in comparison to the subsurface iso-gravity contours of the modeled scenarios. 

25 

In one embodiment, through post processing procedures, the artificially intelligent 
CSS may learn from each successive drilling operation, thus allowing the system to 
develop adaptive behavior routines that will lead to greater accuracy in producing future 
scenarios with each successive drilling operation. The CSS may basically follow a 
30 similar iterative learning process to that of the development of a race car driver from the 
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time he first gets behind the wheel of a car with his driver's education instructor to the 
time he wins the Indianapolis 500. 

In one embodiment, a drill bit is steered based on sampling the gravity associated 
5 with the geology's density variations in real time and automatically making course 
corrections based on a comparison with subsurface iso-gravity contours of geological 
scenario models. It is anticipated that this type of use of GMWD may bring about a 
paradigm shift in directional drilling control. 

10 In other, general embodiments, a drill bit is steered based on any gravity 

measurements made while drilling. Using gravity measurements to steer an apparatus is, 
of course, not limited to drill bits and can include other apparatuses as well. For 
example, real or near-real-time gravity measurements can be utilized to direct any 
subsurface apparatus that requires some type of steering, course correction, or course 

15 verification. 

* * * 

The following examples are included to demonstrate specific embodiments of this 
20 disclosure. It should be appreciated by those having ordinary skill in the art that the 
techniques disclosed in the examples that follow represent techniques discovered by the 
inventor to function well in the practice of the invention, and thus can be considered to 
constitute specific modes for its practice. However, those having ordinary skill in the art 
should, in light of the present disclosure, appreciate that many changes can be made in 
25 the specific embodiments which are disclosed and still obtain a like or similar result 
without departing from the spirit and scope of the invention. 
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Examples 

The following examples emphasize and re-iterate portions of the detailed 
description above, providing a (non-limiting) procedure one may use to accomplish the 
techniques of this disclosure. 

Once the location and proposed subsurface path for a well has been decided upon, 
exemplary embodiments of the present disclosure can involve one or more of the 
following steps. 

1. Produce a best-estimate subsurface three-dimensional density model based on 
geological models developed from geology, geophysics, engineering, and/or well 
data. Such a predicted density model can give the best estimate as well as a range 
of probable subsurface geometries and densities in the form of a probability 
envelope. 

2. Create a three-dimensional grid that will cover the subsurface region of interest, 
including a proposed well path. Throughout the three-dimensional grid, calculate 
the predicted vertical component of gravity due to the predicted density model 
obtained in the previous step. This gravity model can include an error envelope 
that corresponds to the error envelope obtained in the previous step. The 
predicted gravity model that was developed in the previous step can be entered 
into a data processing and modeling computer software, which can be termed a 
Comparative Scenario System (CSS). 

3. A Gravity Measurement While Drilling (GMWD) instrument can be added to a 
drilling unit in operative relation with the drill bit (e.g., above the drill bit or, 
generally, coupled to or near the drill bit). Then, a well is drilled. Each time or at 
selected times during which the drilling effort is stopped (e.g., when it is stopped 
to add another length of drilling pipe), the GMWD unit can measure the vertical 
component of gravity and the dip angle of the well. The drilling effort can also be 
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temporarily stopped whenever an extra set of gravity readings is required or 
desired. A GMWD unit contains one or more gravity sensors designed to operate 
in a drilling environment. Because gravity defines the vertical direction within 
the earth, a gravity meter should be leveled to high precision. The angle between 
vertical as defined by the gravity meter and the axis of the drill pipe is the dip 
angle at the measurement point. 

The gravity data and the dip angle data can be transmitted (e.g., via acoustic 
vibrations of the drilling mud) up the well to the CSS, which may be located at 
the well surface. In addition, the total length of the pipe strung in the well can 
also be entered into the CSS. The dip angle as a function of depth and the total 
pipe length can be used to determine the location of each of the GMWD 
measurements, and hence the location of the drill hole with depth. CSS may be 
an artificial intelligence computer software system that uses a combination of 
gravitational potential forward modeling, and possibly inverse methods, as well as 
numerical analysis techniques known in the art to rapidly evaluate the GMWD 
gravity data by comparing it to the predicted three-dimensional gravitational field 
that was developed in step 2. 

The CSS will make a decision. Using the predetermined tolerances set by the 
predicted gravity error envelope, it can determine if the drill bit is on course or 
needs a course correction. It can also determine if the predicted model is 
sufficiently incorrect to require human intervention. 

If the location of the drill bit is within the pre-determined tolerances, the CSS can 
then calculate a drilling course correction that will automatically be implemented 
through the drilling down-hole guidance system. As long as the GMWD data 
stays within a gravity error envelope, the system may continue automatically. 
The CSS input data, output predicted verses observed correlations and the output 
course corrections can also be transmitted via, e.g., the Internet in real or near real 
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time to anybody that needs to monitor the drilling process, such as the exploration 
team or other decision makers (The Team) that might be located back in an office 
somewhere. Note that the more known and suspected geological and/or target 
information that can be placed into the predicted models, the greater will be the 
5 value of the GMWD data to guide the drilling effort. 

7. Based on the data that is transmitted back to The Team, a human decision could 
be made at any time to intervene and send back, via, e.g., the Internet to the CSS 
at the drill site, a drilling course correction. The CSS may then pass the 

10 correction to a drilling down-hole guidance system for implementation. 

8. If, in step 5, the CSS finds that the predicted gravity scenarios prove to be 
substantially incorrect, the system may inform The Team and ask for a correction 
to the predicted gravity model. The system permits The Team to intervene and 

15 make remotely- accessed, real-time course corrections to the drilling operation 

based on comparative analysis of the actual data obtained near the drill bit in 
comparisons to the predicted gravity models with error envelopes. The CSS can 
assist in obtaining and/or implementing such course corrections. 

20 The CSS may be designed such that, through post processing procedures, the CSS 

will learn from each successive drilling operation, thus allowing the system to develop 
adaptive behavior routines that will lead to greater accuracy in producing future scenarios 
with each successive drilling operation. 
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